the synthesized mass of 56 Ni, M( 56 Ni), to the progenitor metallicity starting 
from Pre-Main Sequence models with masses M = 2 — 7 M and metallicities 
Z = 10~ 5 — 0.10. The interplay between convective mixing and carbon burning 
during the simmering phase eventually rises the neutron excess, 77, and leads 
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1. Introduction 



In addition to the mass, metallicity is one of the few progenitor attributes that can 
leave an imprint on the observational properties of SNIa by affecting the synthesized mass 
of 56 Ni, with important consequences for their use as cosmological standard candles. Up to 
now, attempts to measure Z directly from supernova observat ions have been scarce and their 
results uncertain ( iLentz et al.ll2000l ; iTaubenberger et al.ll2008l ). Measuring Z from the X-ray 
emission of superno ya remnants is a prom ising alternative but as yet has been only applied to 
a single supernova (IBadenes et al.ll2008l ). An al ternative venue is to es timate the supernova 
metallicity as the mean Z of its environment ( IBadenes et al.l 120091 1 . lHamuy et al.l ( 20001 ) 
looked for galacti c age or me t al con tent correlations with SNIa luminosity, but their results 
were ambiguous. Ellis et al.l ( 120081 ) looked for systematic trends of SNIa UV spectra with 
metallicity of the host galaxy, a nd found that t he sp ectral variations were much larger than 
predicted by theoretical models. ICooper et al.l (120091 ) . using data from the Sloan Digital Sky 
Survey and Supern ova Survey concluded t hat prompt SNIa are more luminous in m etal-poor 
systems. Recently, iGallagher et al.l (120081 hereafter G08) and iHowell et al.l (120091 hereafter 
H09), using different methodologies to estimate the metallicity of SNIa hosts, arrived to 
opposite conclusions with respect to the dependence of supernova luminosity on Z. 

There is a long history of numerical simulations of SNIa aimed at predicting the im- 
pact of metallicity and explosive neutronization on their yield s (e.g. 



Bravo et al. 



1992 



Brachwitz et al.ll2000l ; iTravariio et al.l 120051 : IBadenes et al.l 120081 ). iDommguez et al.l ( 12001 
hereafter DHS01) found that the offset in the calibration of supernova magnitudes vs light 
curve (LC) widths is not monotonic in Z and remains smaller than 0.07 m for Z < 0.02. 
Kasen et al.l (120091 ) concluded that the width-luminosity relation ship depends weakly on the 
metallicity of the progenitor. From an analytical point of view, iTimmes et al.l ( 120031 . here- 
after TBT03) using arguments from basic nuclear physics predicted a linear relationship 
between M( 56 Ni) and Z. The conclusions of TBT03 relied on two main assumptions: first, 
that most of the 56 Ni is synthesized in material that burns fully to NSE and, second, that 
a fiducial SNIa produces a mass M^j ~ .6 Mp of Fe-group nuclei whose p is not modi- 
fied during the explosion. iPiro fc Bildstenl (120081 ) and IChamulak et al.l ( 120081 ). based on the 
same assumptions as TBT03, extended their analysis taking into account the neutronization 
produced during the simmering phase. 



In this paper, we show that the first assumption of TBT03 does not hold for most 
SNIa. Indeed, for a SNIa that produces Mp° ~ 0.6 M the fractio n of 56 Ni synthesized 
out of NSE exceeds ~ 30%. With respect to the second assumption, iMazzali et al.l ( 120071 . 
hereafter M07) showed, based on observational results, that the mass of Fe-group nuclei 
ejected by SNIa spans the range from 0.4 to 1.1 M . This range cannot be accounted for by 
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metallicity variation within reasonable values. Accordingly, our working hypothesis is that 
the yield of 56 Ni in SNIa is governed by a primary parameter different from Z . In our one- 
dimensional models the primary parameter is the DDT density, Pddt ; although in nature 
it may be something else such as the expansion rate during the deflagration phase. The 
initial metallicity is a secondary factor that can give rise to scatter in the value of M( 56 Ni) 
either directly (linear scenario), by affecting the chemical composition of the ejecta for a 
given value of the primary parameter, or indirectly (non-linear scenario), by modifying the 
primary parameter itself. The understanding of which one of these two characters is actually 
being played by Z is of paramount importance. 



2. The effect of metallicity on the yield of Ni 



We explore the sensitivity of M( 56 Ni) to the progenitor metallicity, starting from Pre- 
Main Sequence models of masses, M , in the range 2 — 7 M Q and metallicities, Z, from 10~ 5 
to 0.10, as given in the first column of Table [TJ The initial m ass fracti o ns of all the isotopes 
with A > 6 have been fixed in solar proportion, according to iLodderd (120031 ); consequently, 
we adopt for the solar metallicity the value Z Q = 0.014. Each presupernova model has been 
evolved from the Pre-Main Sequence to the Thermal Pulse (TP) AGB phase, in order to 
determine the mass, M core , and chemical structure of the C-0 core left behind. Afterwards, 
an envelope of the appropriate size to reach the Chandrasekhar mass, Mch, has been added 
on top of the C-0 cores, and these structures have been fed as initial models to a supernova 
hydrocode. Finally, the explosive nucleosynthesis has been obtained with a post-processing 
nucleosynthetic code. 



The hydrostatic evolution has been computed by means of the FRANEC code (jChieffi et al 



19981 ). With respect to the calculations of DHS01, the code has been updated in the input 
physics. For the purposes of the present paper, the most i mportant changes c oncern the 
12 C (a, 7) 16 O reaction rate, which is calculated according to iKunz et al.l ( 2002 ) instead of 
Caughlan et al. ( 1985h . and the treatm ent of convective mixing during the late part of the 



core- He burning ( 



Straniero et al.ll2003l ). 



The presupernova model is a Chandrasekhar mass WD built in hydrostatic equilibrium 
with a central density p c = 3 x 10 9 g cm -3 . The composition of the envelope of mass 
Mch — M COTe is the same as that of the outermost shell of the C-0 core. Thus, instead 
of assuming C/0=1, as in DHS01, we adopt the C/O ratio obtained as a result of He- 
shell burning during the AGB phase. The effect of changing p c and the composition of the 
envelope has been tested in several models, as explain ed later. We leave aside other eventual 
complexities of pre-supernova physics like rotation (IPiersanti et al.l 120031 ; lYoon fc Langer 
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2004h . 



The internal composition of the WD is eventually modified during the simmering phase, 
due to the combined effects of convective mixing, carbon burning and electron captures. The 
first two phenomena affect the carbon abundance within the core, while the latter leads to an 
increase of rj. The average (within the WD) carbon consumption and neutron excess increase 
during the simmering phase are A"K( 12 C) ~ —1.66 x 10~ 3 mol g~ x and An = — |AF( 12 C) 
( IChamulak et al.ll2008l ). We assume that convective mixing is limited to the C-0 core, which 



implies that the change in the neutron excess within the core is An m 1.11 x 10~ 3 Mch/M core . 
We have also exploded several models disregarding the simmering phase, to which we will 
refer in the following as stratified models. 



The supernova hydrodynamics code we have used is the same as in lBadenes et all (12003 ). 
As in DHS01, the present models are based on the delayed-detonation paradigm (IKhokhlov 
19911 ). To address the linear scenario we take pddt independent of Z. In this case, pddt = 
3 x 10 7 g cm -3 , although simulations with pddt in the range 1 
reported. 



3 x 10 7 g cm 3 are also 



For the non-linear scenario we have adopted the criterion that a DDT is induced when 
the laminar flame thickne ss, <5i am , becomes of the order of the turbulent Gibson length Iq 
(IRopke fc Niemeyerl 120071 ). w ith the flame p roperties ( velocity and width) de pending on the 
abundan ces of 12 C (eq. 2 2 in IWooslevi 120071 ) and 22 Ne (IChamulak et al.ll2007l ). and hence on 
Z and rj. iTownsley et al.l ( 120091 ) concluded from 2D simulations of SNIa that the metallicity 
does not affect the dynamics of the explosion, and so the turbulence intensity is independent 
of Z. Thus, for a given turbulent intensity a change in Z can be compens ated by a change 
in pn nT in order to recover the condition 5\ am / Iq ~ 1 (see the discussion in IChamulak et al. 
20071 ). In this scenario we have scaled the transition density as a function of the local 
chemical composition as follows: 



PDDT ocX( 12 C)- L3 (l 



129?7) 



-0.6 



In order to introduce an rj dependence in the above expression we have assumed, for simplic- 
ity, that the bulk of neutronized isotopes synthesized during the simmering phase accelerates 
the carbon consumption rate the same way 22 Ne does. 



2.1. Presupernova evolution 

The results of the hydrostatic evolution of our presupernova models are shown in Tabled] 
For each Mq and Z we give: M core , the central abundance of 12 C and rj in stratified models, 
X C ( 12 C) and rj c , and the same quantities in the models accounting for the simmering phase, 
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X s i m ( C) and 77 s i m . In comparison with DHS01, the present models span a larger range of 
Z, as DHS01 computed models with Z < Z & . In the common range of Z and M the results 



are comparable, although the adopted rate of the C (a, 7) O reaction leads to a slightly 
larger carbon abundance than in DHS01. The differences in M core between our models and 
those of DHS01 are smaller than 0.06 M Q . The central carbon to o xygen ratio and M C „ TC . we 
obtain, and their dependencies with Z and Mq, agree as well with lUmeda et al.l (119991 ) . 



2.2. Mass of M( 56 Ni) ejected 

The results of the explosion simulations are summarized in Figs. [T] and Figured] shows 
the dependence of M( 56 Ni) on Z. For the stratified models, we obtain the same range of 
variation of M( 56 Ni) with respect to M at given Z as DHS01: 0.06 M , although the 56 Ni 
yields do not match with DHS01 because they used different values of p c = 2 x 10 9 g cm -3 
and Pddt = 2.3 x 10 7 g cm -3 . The models accounting for the simmering phase behave like 
the stratified models with respect to variations in M Q and Z, although with a smaller total 
M ( 56 Ni) due to electron captures during the simmering phase. The dependence of M( 56 Ni) 
on Z can be approximated by a linear function: 

M( 56 Ni) oc f(Z) = 1 - 0.075 J- , (2) 

Zq 

while stratified models can be approximated by: M( 56 Ni) oc 1 — 0.0Q9Z/Z & , i.e. the slope 
of the linear function is quite insensitive to the carbon simmering phase. 

To explore the non-linear scenario we have computed models accounting for the sim- 
mering phase with fixed M = 5 M . Introducing a composition dependent pddt produces 
a qualitatively different result because the relationship between M( 56 Ni) and Z is no longer 
linear, especially at low metallicites for which a larger Pddt is obtained, implying a much 
larger M( 56 Ni). Our results can be fit by a polinomial law: 

M( 56 Ni) oc f(Z) = 1 - 0.18-£ (l - °- 10 ^) • ( 3 ) 

Both the central density at the onset of thermal runaway and the final C/O ratio in the 
accreted layers have a minor effect on M( 56 Ni) within the explored range. 

TBT03 proposed a linear relationship between M( 56 Ni) and Z: M( 56 Ni) oc 1-0.057 Z/Z Q 
(dotted line in Fig. [T]). In all of our present models we find a steeper slope. The reason for 
this discrepancy lies in the assumption by TBT03 that most of the 56 Ni is synthesized in 
NSE. In our models a sizeable fraction of 56 Ni is always synthesized during incomplete Si- 
burning, whose final composition has a stronger dependence on Z than NSE matter. As 
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Hix fc Thielemannl ( 119961 ) showed, the mean neutronization of Fe-peak isotopes during in- 



complete Si-burning is much larger than the global neutronization of matter because neutron- 
rich isotopes within the Si-group are quickly photodissociated, providing free neutrons that 
are efficiently captured by nuclei in the Fe-peak group, favouring their neutron-rich isotopes. 
Figure [2] shows that up to ~ 60% of Mp° can be made out of NSE, the actual fraction de- 
pending essentially on the total mass of Fe-group elements ejected. Thus, the less M( 56 Ni) 
is synthesized, the larger fraction of it is built during incomplete Si-burning and the stronger 
is its dependence on Z. 



3. Discussion 

The results presented in the previous section show that the metallicity is not the pri- 
mary parameter that allows to reproduce the whole observational scatter of M( 56 Ni), for 
a reasonable range of Z. We have also shown that a possible dependence of the primary 
parameter on Z, would lead to a non-linear relationship between M( 56 Ni) and Z, as in Eq. El 
However, as we will show in the following, it would be possible to unravel the way M( 56 Ni) 
depends on Z by means of future accurate measurements of SNIa properties. 

We start analysing the amount of the scatter induced by the dependence of M( 56 Ni) on 
Z given by Eq. [3J For simplicity we follow the procedure of M07 to estimate the supernova 
luminosity and LC width. The peak bolometric luminosity, L, is determined directly by the 
mass of 56 Ni synthesized (in the following, all masses are in M and energies are in 10 51 ergs): 

L [M( 56 Ni)j = 2 x 10 43 M( 56 Ni) erg s" 1 , (4) 

while the bolometric LC width, r, is determined by the kinetic energy, E^, and the opacity, 
k: t (x K 1/ ' 2 i? k r 1//4 . The kinetic energy is given by the difference of the WD initial binding 
energy, |BE|, and the nuclear energy released, the latter being related to the final chemical 
composition of the ejecta: E k » 1.56M( 56 Ni) + 1.74 [M Fe - M( 56 Ni)] + 1.24Mi ME - |BE| 
, where M Fe is the total mass of Fe-group nuclei and Mime is the mass of intermediate- 
mass elements (IME). The opacity is provided mainly by Fe-group nuclei and IMEs: k oc 
M Fe + O.IMjme- We have taken |BE| = 0.46, which is a good approximation given the small 
variation of binding energy with initial central density: |BE| is in the range 0.44 — 0.47 for 
p c = 2 — 4 x 10 9 g cm -3 . To reduce the number of free parameters we further link Mime 
to M Fe imposing that the ejected mass is the Chandrasekhar mass (Mch ~ 1-38 M Q in our 
models), and that the amount of unburned C+O scales as Mqo ~ 0-3Mj 2 ME , as deduced from 
our models. Thus, M Fe + Mime + 0.3Mj 2 ME = M Ch . Furthermore, the mass of 56 Ni is linked 
to the mass of Fe-group nuclei by M( 56 Ni) = M£° x f(Z) = (M Fc - M cc ) x f(Z), where 
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f(Z) is given by Eq. [3] or a similar function, and M ec is the mass of the neutron-rich Fe- 
group core (due to electron captures during the explosion). We have taken M ec « 0.14 M , 
which is representative of the range of masses obtained in our models: 0.10 — 0.16 M Q for 
p c = 2 — 4 x 10 9 g cm" 3 . Finally, to compare with observed values a scale factor of 24.4 is 
applied to the value of r thus obtained, as in M07. Putting all these together, we obtain the 
following expression for the bolometric LC width (in days) function of M( 56 Ni) and Z: 



1/2 



r [M( 56 Ni),Z] = 21.9' 



1/4 



| (j^ ~ 0.115) M(»Ni) - 1.46 + 2.09^1 - 0.482^^1 



(5) 

The relationship between L and r derived from Eqs. [3], H] and [5] is displayed in Fig. [3] for 
three different metallicities along with observational data. There are also represented the 
relationships obtained by substituting Eq. [3] by the M( 56 Ni) vs Z dependences proposed 
by TBT03 and Eq. 5 in H09. Our Eq. [3] gives a wider range of M( 56 Ni), which accounts 
better for the scatter of the observational data. Indeed, if real SNIa follow Eq. [3], deriving 
supernova luminosities from Z- uncorrected LC shapes might lead to systematic errors of up 
to 0.5 magnitudes. 

To estimate the bearing that the metallicity dependence of M( 56 Ni) can have on cos- 
mological studies that use a large observational sample of supernovae, we have generated a 
virtual population of 200 SNIa that has been analyzed following the same methodology as 
G08 and H09. Each virtual supernova has been randomly assigned a progenitor metallicity, 
from a uniform distribution of log(Z) between Z m[n = 0.1Z Q and Z max = 3Z Q , and an M Fe , 
uniformly distributed in the range from 0.31 to 1.15 M . The minimum and maximum 
M( 56 Ni) thus obtained (computed with Eq. [3] and M ec = 0.14 M ) are 0.1 and 1 M Q , and 
the bolometric LC width, r, lies in the range 15 — 24 days. A Z-uncorrected mass of 56 Ni, 
M^, has then been obtained as the value of M( 56 Ni) that would give the same r if Z = Z & . 
The Mf & so computed gives an idea of the effect of fitting an observed SNIa LC with a 
template that takes no account of the supernova metallicity. From Eq. HI we estimate the 
Hubble Residual, HR, of each virtual SNIa at: HR = 2.5 log (M|/M( 56 Ni)). As a final step 
we have added gaussian noise with a = 0.1 to both HR and \og(Z), to simulate the effect of 
observational uncertainties. 

A linear relationship HR = a + (3log(Z) has then been fit to the noisy virtual data 
by the least-squares technique, as in G08 and H09. Figure 0] shows the results for 10,000 
realizations of the noisy virtual dataset. The histogram gives the number counts of the slope 
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in the 10,000 realizations. The whole process has been repeated by using Eq. [2] (i.e. the 
linear scenario) to represent the dependence of M( 56 Ni) on Z and the results are also shown 
in Fig. |H From the Figure it is clear that, for a large enough set of SNIa whose luminosity 
and metallicity are measured with small enough errors, it is possible to discriminate between 
the linear and non-linear scenarios. In our numerical experiment, the mean value of (3 is 0.13 
in the first case and 0.26 in the second case, both with a standard deviationof 0.02. 

Figure H] shows also the observational results obtained by G08, who approximated the 
metallicities of the SNIa in their sample by the Z of the host galaxy, obtained from an 
empirical galactic mass-metallicity relationship. The striking match between our results 
based on the non-linear scenario and those of G08 must be viewed with caution in view 
of the observational uncertainties involved in measuring supernova metallicities and the 
limitations of our models (i.e. the assumption of spherical symmetry). Recently, using a 
different method of determination of the SNIa metallicity, H09 arrived to a result opposite 
to that of G08, i.e. they found that HR is uncorrelated with Z, leading to a distribution 
centered around j3 ~ 0. Thus, until such discrepancies are resolved it is not possible to 
draw any firm conclusion about the metallicity effect on SNIa luminosity. However, it is 
worth stressing that the simultaneous measurement of supernova luminosity and metallicity 
for a large SNIa set would strongly constrain the physics of the deflagration-to-detonation 
transition in thermonuclear supernovae, one of the key standing problems in supernova 
theory. 

This work has been partially supported by the MEC grants AYA2007-66256 and AYA2008- 
04211-C02-02, by the European Union FEDER funds, by the Generalitat de Catalunya, and 
by the ASI-INAF 1/016/07/0. CB thanks Benoziyo Center for Astrophysics for support 

REFERENCES 

Badenes, C, Bravo, E., Borkowski, K. J., & Dommguez, I. 2003, ApJ, 593, 358 

Badenes, C, Bravo, E., & Hughes, J. P. 2008, ApJ, 680, L33 

Badenes, C, Harris, J., Zaritsky, D., & Prieto, J. L. 2009, ApJ, 700, 727 

Brachwitz, F., Dean, D. J., Hix, W. R., et al. 2000, ApJ, 536, 934 

Bravo, E., Isern, J., Canal, R., & Labay, J. 1992, A&A, 257, 534 

Caughlan, G. R., Fowler, W. A., Harris, M. J., & Zimmerman, B. A. 1985, Atomic Data 
and Nuclear Data Tables, 32, 197 



-9 - 

Chamulak, D. A., Brown, E. F., & Timmes, F. X. 2007, ApJ, 655, L93 

Chamulak, D. A., Brown, E. F., Timmes, F. X., & Dupczak, K. 2008, ApJ, 677, 160 

Chieffi, A., Limongi, M., & Straniero, O. 1998, ApJ, 502, 737 

Contardo, G., Leibundgut, B., & Vacca, W. D. 2000, A&A, 359, 876 

Cooper, M. C, Newman, J. A., & Yan, R. 2009, ApJ, 704, 687 

Dommguez, L, Hoflich, P., & Straniero, O. 2001, ApJ, 557, 279 

Ellis, R. S., Sullivan, M., Nugent, P. E., et al. 2008, ApJ, 674, 51 

Gallagher, J. S., Garnavich, P. M., Caldwell, N., et al. 2008, ApJ, 685, 752 

Hamuy, M., Trager, S. C, Pinto, P. A., et al. 2000, AJ, 120, 1479 

Hix, W. R. & Thielemann, F. 1996, ApJ, 460, 869 

Howell, D. A., Sullivan, M., Brown, E. F., et al. 2009, ApJ, 691, 661 

Kasen, D., Rdpke, F. K., & Woosley, S. E. 2009, Nature, 460, 869 

Khokhlov, A. M. 1991, A&A, 245, 114 

Kunz, R., Fey, M., Jaeger, M., et al. 2002, ApJ, 567, 643 

Lentz, E. J., Baron, E., Branch, D., Hauschildt, P. H., & Nugent, P. E. 2000, ApJ, 530, 966 
Lodders, K. 2003, ApJ, 591, 1220 

Mazzali, P. A., Rdpke, F. K., Benetti, S., & Hillebrandt, W. 2007, Science, 315, 825 

Phillips, M. M., Krisciunas, K., Suntzeff, N. B., et al. 2006, AJ, 131, 2615 

Piersanti, L., Gagliardi, S., Iben, I. J., & Tornambe, A. 2003, ApJ, 598, 1229 

Piro, A. L. & Bildsten, L. 2008, ApJ, 673, 1009 

Rdpke, F. K. & Niemeyer, J. C. 2007, A&A, 464, 683 

Stanishev, V., Goobar, A., Benetti, S., et al. 2007, A&A, 469, 645 

Straniero, O., Dommguez, I., Imbriani, G., & Piersanti, L. 2003, ApJ, 583, 878 

Taubenberger, S., Hachinger, S., Pignata, C, et al. 2008, MNRAS, 385, 75 



- 10 - 

Timmes, F. X., Brown, E. F., & Truran, J. W. 2003, ApJ, 590, L83 
Townsley, D. M., Jackson, A. P., Calder, A. C, et al. 2009, ApJ, 701, 1582 
Travaglio, C, Hillebrandt, W., & Reinecke, M. 2005, A&A, 443, 1007 
Umeda, H., Nomoto, K., Yamaoka, H., & Wanajo, S. 1999, ApJ, 513, 861 
Wang, X., Li, W., Filippenko, A. V., et al. 2009, ApJ, 697, 380 
Woosley, S. E. 2007, ApJ, 668, 1109 
Yoon, S. & Langer, N. 2004, A&A, 419, 623 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 11 - 



Table 1. Properties of CO cores at the beginning of TPs 





M 
(M©) 


-^core 

(M©) 


X c ( 12 C) 


Vc 




^sim ( 12 C) 


^7sim 




10~ 5 (0.23) 


3 


0.801 


0.28 


8.7 x 10" 


-7 


0.38 


1.9 x 10' 


-3 


10~ 5 (0.23) 


5 


0.903 


0.25 


8.6 x 10- 


-V 


0.31 


1.7 x 10' 


-3 


10~ 5 (0.23) 


6.5 


1.052 


0.20 


8.4 x 10' 


-7 


0.24 


1.5 x 10' 


-3 


0.014 (0.269) 


3 


0.615 


0.23 


1.2 x 10- 


-3 


0.36 


3.7 x 10" 


-3 


0.014 (0.269) 


5 


0.848 


0.32 


1.2 x 10- 


-3 


0.38 


3.0 x 10' 


-3 


0.014 (0.269) 


7 


1.005 


0.22 


1.2 x 10- 


-3 


0.26 


2.7 x 10- 


-3 


0.040 (0.31) 


3 


0.608 


0.18 


3.4 x 10- 


-3 


0.33 


6.0 x 10- 


-3 


0.040 (0.31) 


5 


0.843 


0.32 


3.4 x 10- 


-3 


0.39 


5.4 x 10- 


-3 


0.040 (0.31) 


7 


1.032 


0.28 


3.4 x 10- 


-3 


0.30 


5.1 x 10- 


-3 


0.10 (0.38) 


3 


0.624 


0.15 


7.3 x 10- 


-3 


0.29 


10.1 x 10 


-3 


0.10 (0.38) 


5 


0.846 


0.29 


7.6 x 10" 


-3 


0.34 


9.6 x 10" 


-3 


0.10 (0.38) 


7 


0.962 


0.23 


7.6 x 10- 


-3 


0.26 


9.5 x 10- 


-3 



a Initial metallicity (and helium abundance) 
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Ni yield vs initial metallicity for different initial masses and explosion parameters. 



Black triangles: stratified models with p c = 3 x 10 9 g cm -3 and Pddt = 3 x 10 7 g cm -3 . Red 
(filled) circles: models exploded with the same parameters but accounting for the simmering 
phase. Blue (empty) circles: the same as red circles except that pddt is a function of the 
local (at flame position) Z through X( 12 C) and r\ (Eq. [T|) . Blue crosses and blue stars: 
the same as blue circles except that p c = 2 x 10 9 g cm -3 (crosses) or p c = 4 x 10 9 g cm -3 
(stars). Blue squares: the same as blue circles except that the composition of the envelope 
is composed by equal amounts of carbon and oxygen, i.e. C/O = 1 as compared to values 
ranging from C/O = 1.5 to 2.3 as taken from the He-shell burning during the TP phase. The 
lines represent the linear relationship between M( 56 Ni) and Z proposed in TBT03 (dotted 
line), and our Eq. [3] (dashed line). 
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Fig. 2. — Fraction of Fe-group elements (from Ti to Zn) synthesized in layers experiencing 
incomplete Si-burning (T max < 5.2 x 10 9 K) as a function of the total mass of Fe-group 
elements ejected. The filled square is a reference model with Pddt = 3 x 10 7 g cm -3 , 
Z = 10~ 5 , and Mo = 5M , that accounts for the simmering phase. The empty squares 
belong to models with varying pddt, from 2.6 x 10 7 to 10 7 g cm -3 . All the above models 
are linked by a solid line to help guiding the eye. The rest of models show the sensitivity to 
different model parameters with respect to the reference model, as follows: crosses, varying 
Z; empty circles, varying M ; stars, varying p c 
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Fig. 3. — Sensitivity of the bolometric light curve luminosity-width relationship to the initial 
metallicity. The solid lines belong to Eq. [3]of the present work for Z = 0.1 Z & (green), Z = Z Q 
(blue), and Z = 3Z & (red). There are also shown the curves obtained using the M( 56 Ni) vs 
Z rel ationships of TBT03 ( dotted lines) and H09 (dashed lines). The stars represent SNIa 
from IContardo et al.l (120001 ) . whose light-curve width h as been com puted as in M0 7, whil e 
the asterisks re present SN1999ac Jphillips et al.l l20Q6h . SN2003du Jstanishev etldl l2007h . 
and SN2005df (IWang et al.l 120091 ) (by order of decreasing LC width), estimated from the 
published bolometric light curves 



- 15 - 




Fig. 4. — Number distribution of the slope, /?, of the least-squares fit to the Hubble residual 
vs metallicity. The statistics used 10,000 data sets, each one obtained adding gaussian noise 
(<r = 0.1 in HR, and 0.1 dex in Z) to a virtual random population of 200 SNIa generated 
using either the quadratic M( 56 Ni) vs Z relationship (Eq. |3l horizontally hatched histogram) 
or the linear relationship (Eq. |2j vertically hatched histogram). The vertical dashed line 
shows the slope measured by G08 



